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Abstract 

A nonperturbative correction to the thermal nucleation rate of critical bubbles in a 
first order phase transition is estimated. The correction originates from large-amplitude 
fluctuations which may be present before the transition occurs. Using a simple model 
of a scalar field in a double-well potential, we present a method to obtain a corrected 
potential which incorporates the free-energy density available from large-amplitude 
fluctuations, which is not included in the usual perturbative calculation. For weaker 
phase transitions, the nucleation rate can be much larger than the rate calculated via 
perturbation theory. As an application of our method, we show how nonperturba- 
tive corrections can both qualitatively and quantitatively explain anomalously high 
nucleation rates observed in 2-d numerical simulations. 
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Although the simplest first order phase transitions are characterized by a discontinuous 
jump of a scalar order parameter between two distinct phases, they do not all proceed 
in the same way For very strong first order phase transitions, where the free-energy 
barrier between the phases is large, the transition is initiated by the nucleation of critical- 
sized bubbles of the new phase in the background of the metastable {e.g., super-cooled) old 
phase. By definition, these critical bubbles are just large enough to overcome their surface 
tension and grow, eventually converting the whole medium to the new phase. The large 
barrier between the two phases suppresses large- amplitude thermal fluctuations of the order 
parameter; an initial metastable state is well-defined, as no fraction of the volume is in the 
new phase before the transition occurs. In this case, the metastable phase can be regarded 
as "homogeneous", as only very small-amplitude thermal fluctuations are present. This is 
the situation described by hanger's theory of homogeneous nucleation 0, or, in the context 
of relativistic quantum field theories, by the work of Coleman and Callan 0. 

Besides the decay of the "near-homogeneous" metastable state described by nucleation 
theory, one can investigate the evolution of an unstable initial state which is characterized by 
considerable phase mixing. Within the context of condensed matter systems, this situation 
corresponds to a quench within the unstable "spinodal" region of the two-phase diagram. In 
this case, the two phases separate by the mechanism known as "spinodal decomposition"; 
small- amplitude, long-wavelength fluctuations grow exponentially fast, forming domains of 
the two phases which will eventually coarsen, as the system approaches its final equilibrium 
state. 

In this letter we will address the dynamics of phase transitions characterized by an 
initial state which lies within the "grey zone" between homogeneous nucleation and spinodal 
decomposition. Looking at the whole "spectrum" of first order phase transitions, from 
very strong to very weak, it is clear that the amount of phase-mixing of the initial state will 
strongly influence the subsequent dynamics of the transition. However, the standard method 
of calculating the nucleation rate employs Gaussian perturbation theory, which is valid only 
for small amplitude fluctuations 0]. For strong transitions this approximation is valid. But 
for weaker transitions, large amplitude fluctuations are more abundant, and can have an 
important effect. Our goal is to present an approximate method by which the presence of 
large-amplitude fluctuations is consistently incorporated into the calculation of nucleation 
rates. Thus, we are implicitly assuming that we are close enough to the regime described 
by homogeneous nucleation that we can still distinguish between the two low-temperature 
phases. 

Large-amplitude thermal fluctuations will be modelled by the so-called sub-critical bubble 
method . Recent results have shown that modelling the dominant fluctuations by sub- 
critical bubbles is in excellent agreement with 3-d simulations The model utilizes the 
fact that along with the nucleation of critical bubbles in the meta-stable phase, smaller 
size, though still large amplitude, "sub-critical" bubbles will also be nucleated (and in much 
greater number because they have a lower free energy). These bubbles by definition will 
always shrink and eventually disappear, but there will always be some non-zero equilibrium 
number density Usb at a given temperature. Their presence may lead to large corrections on 
nucleation rates. 

To begin, let us consider the standard model of a phase transition, in which the order 
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parameter is a real scalar field 0, which has a quartic double-well potential of the form 

^(0) = \rn^' - + A04_ (1) 

This potential has two minima, one at = and at = 0+, which represent the two phases 
of the system. It can be thought of as the homogeneous part of a typical phenomenological 
Ginzburg-Landau coarse-grained free-energy density (the cubic term can always be made into 
a linear term), or as some effective potential where additional degrees of freedom coupled to 
have been integrated out. Our analysis will be purely classical, valid for T ^ m, where 
m is the mass of the low-energy mesonic excitations in the associated quantum theory. All 
relevant field configurations contain many quanta. 

We would like to incorporate the free-energy density associated with large-amplitude, 
nonperturbative fluctuations into the computation of the decay rate. In the spirit of the 
renormalization group approach, this should be equivalent to an effective "coarse-graining" 
of the classical potential; averaging over these large-amplitude fluctuations will lead to a shift 
in the background free-energy density and decay barrier, which in principle can be translated 
into a change in the bare couplings of the model. We can understand how to estimate the 
effective coarse-graining by first studying the thin-wall limit of critical bubble nucleation. 

In the standard theory, which neglects phase mixing, the nucleation rate F is proportional 
to e~^'=^^'^, where F^h is free energy needed to form a critical bubble in the metastable back- 
ground. For an arbitrary thin- walled spherical bubble of radius R and amplitude 0thin ^ 0+, 
where thin-walled means the radius R is much greater than the bubble wall thickness, the 
free energy of the bubble takes the well-known form p] 

F,UR) = A7rR^a-—R'AV. (2) 

This formula has a simple physical interpretation. The first term is the energy it costs to 
form the bubble wall, where a = | / dr{d(f)/dr)'^ is the surface tension. The second term is 
the energy "gained" by converting a spherical volume of the metastable phase into the lower 
energy phase. Therefore, is defined as the difference in free-energy density between 
the background medium and the bubble's interior. Since 0thin ^ 0+, for a homogeneous 
background (metastable) we can write, 

AK, = nO)-V^(0+), (3) 

where we have explicitly used the subscript to stress that this is for the case with no phase 
mixing. 

If there is significant phase mixing in the background metastable state, its free-energy 
density is no longer V^(0). One must also account for the free-energy density of the non- 
perturbative, large-amplitude fluctuations. Since there is no formal way of deriving this 
contribution outside improved perturbative schemes, we propose to estimate the corrections 
to the background free-energy density by following another route. We start by writing 

free — energy density of metastable state = V{0) + J-'sc, (4) 



3 



where J^sc is the nonperturbative contribution to the free-energy density due to the large 
amphtude fluctuations, which we assume can be modelled by subcritical bubbles. We will 
calculate T^c further below. 

We thus define the effective free-energy difference AV^g, which includes corrections due 
to phase mixing, as 

AKg = Ayo + ^sc (5) 

which is the sum of the free-energy difference calculated in the standard way [eq. (j^)], and 
the "extra" free-energy density due to the presence of subcritical bubbles. Henceforth, the 
subscript 'eg' will stand for "coarse-grained" . 

We note that while we have made a correction to AV, we have not made any correction 
to the surface tension a. Since we are considering the thin-wall limit, as long as (0) is small, 
which is true if subcritical bubbles do not occupy a large fraction of space, the correction 
to (J will be subdominant. [Note that the presence of subcritical bubbles may shift (0) by 
roughly l]ie~^'/^0i, where 0i is the amplitude of a given fluctuation, and Fi its associated 
free energy.] Thus, the arguments here give a lower bound on the magnitude of the correc- 
tions. Later on, both volume and surface corrections will be automatically included in the 
calculation. 

Since a critical size bubble is defined as the bubble for which all forces on the bubble 
wall cancel, i.e. dF/dR\R^^ = 0, we can now use eq. (|^) to obtain the free energy needed to 
form a thin-wall critical bubble in a background with subcritical bubbles 

F., = ^-^RliAVo + ^sc) = (6) 
3 3 (AVo + J^scV 

and the radius of the critical bubble is 

Equations (|^) and (0) warrant several comments. First, in the limit of a very strong phase 
transition, subcritical bubbles are suppressed {J^sc 0), and both Feb and Rch approach the 
standard expressions for the free energy and radius of a critical bubble in a homogeneous 
background. Second, notice that J-'sc acts in the same way as the free-energy difference AVq. 
The presence of subcritical bubbles is equivalent to extra free energy in the medium, which 
enhances the nucleation of critical bubbles. In particular, for potentials near degeneracy such 
that AVo ^ -^sc, the nucleation rate of critical bubbles T ~ e~^'=^^'^, can be much greater 
than in the case ignoring the presence of subcritical bubbles. 

Finally, notice that as AVq — > 0, neither the critical-bubble energy nor its radius become 
infinite. For temperature-dependent potentials which (ignoring the corrections from subcrit- 
ical bubbles) are degenerate at the critical temperature Tc, the nucleation rate F ~ is 
finite. In fact, the nucleation rate of critical bubbles may be non-zero even above the critical 
temperature (again, using the uncorrected expression for the potential). This is a testable 
prediction of our method which, of course, is sensitive to the equilibrium number-density of 
subcritical bubbles. 
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This final comment suggests an important point. Since for degenerate potentials (tem- 
perature dependent or not) no critical bubbles should be nucleated, taking into account 
subcritical bubbles must lead to a change in the coarse-grained free-energy density (or po- 
tential) describing the transition. Thus, it should be possible to translate the "extra" free 
energy available in the system due to the presence of subcritical bubbles in the background 
into a corrected potential for the scalar order parameter. We will write this corrected po- 
tential as Kg(0)- 

The standard coarse-grained free energy is calculated by integrating out the short wave- 
length modes (usually up to the correlation length) from the partition function of the system, 
and is approximated by the familiar form [^] 

Fc, = /c?'rQ(V0)2 + Kg(0)) . (8) 

How do we estimate V"cg? One way is to simply constrain it to be consistent with the thin 
wall limit. That is, as Kg(0) approaches degeneracy {i.e. AVcg{(p) — > 0), it must obey the 
thin wall limit of eq. (^. Note that with a simple rescaling, the potential of eq. (|I|) can be 
written in terms of one free parameter. Thus, the thin wall constraint can be used to express 
the corrected value of this parameter in terms of J^^c in appropriate units. The free energy of 
the critical bubble is then obtained by finding the bounce solution to the equation of motion 
— dVcg{(t>) / d(j) = by the usual shooting method, and substituting this solution into 
eq. d). 

Therefore, in order to determine V^g, we must first calculate the free-energy density J-'sc 
of the subcritical bubbles. As a first step, we follow the work of Ref. [0, to obtain the 
equilibrium number density Ush of subcritical bubbles. If we define the distribution function 
/ = d'^nsh/dRd(j)A, then f{R,(f)A,t)dRd(f)A is the number density of bubbles with radius 
between R and R + dR and amplitude between 0^ and (pA + d(j)A at time t. It satisfies the 
Boltzmann equation, 

dfiR,<PA,t) \ \df 

di " -11^1^ + (1-7)^0^+ 

-/VG-Thcrm - iG+^o- (9) 

The first term on the RHS is the shrinking term (note that v = dR/ dt is negative), the second 
term is the nucleation term where G is the nucleation distribution function, which is defined 
by r = / dRd(f)G, and ro^+ is the nucleation rate per unit volume of subcritical bubbles 
from the "0" phase (the initial phase) to the phase. The division of the system into two 
phases depends on the particular application at hand, as will be clear in the example below. 
By the Gibb's distribution, 6*0-*+ = Ae~^''^^^''^'^^^'^ , where A is a constant independent of R 
and (f). 

The factor 7 is defined as the fraction of volume in the phase, and is obtained by 
summing over subcritical bubbles of all amplitudes within this phase. The third term is a 
phenomenological thermal destruction term (see work by Gelmini and Gleiser in Ref. |^), 
where V is the volume of a bubble of radius R, and GThcrm = o,T/V, where a is a constant. 
The fourth term is the inverse nucleation term. For more details about this Boltzmann 
equation, see Ref. 0, which has improved upon the work of Gelmini and Gleiser (Ref. P). 
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The free energy of the subcritical bubbles is determined by modelhng them as Gaussian 
fluctuations with amplitude (j)A and radius R, 

0sc(r)=0Ae-^'/^'. (10) 

The free energy of a given configuration can then be found by using the general formula, 

i^sb = / d'r (i(V0sc)' + l^(0sc)) . (11) 

Although this approach only includes one particular shape out of all possible field configura- 
tions, the agreement between theory and numerical experiments indicates that the Gaussian 
profile is an adequate ansatz for the dominant large-amplitude thermal fluctuations. 

The equilibrium number density of subcritical bubbles is found by solving eq. (^ with 
df/dt = 0, imposing the physical boundary condition /(r — > oo) = 0. Once we know the 
distribution function and free energy for a bubble of a given radius R and amplitude (pA, we 
can estimate the total energy density of the Gaussian subcritical bubbles, summed over all 
relevant radii and amplitudes. We can write, in general, 

^sc ~ / / ^shTT^K—dRd^A, (12) 

where 0min defines the lowest amplitude within the "+" phase, typically (but not necessarily) 
taken to be the maximum of the double-well potential. -Rmin is the smallest radius for the 
subcritical bubbles, compatible with the coarse-graining scale. For example, it can be a 
lattice cut-off in numerical simulations, or the mean-field correlation length in continuum 
models. As for -Rmax, if is natural to choose it to be the critical bubble radius. 

As an application of the above method, we will investigate nucleation rates in the context 
of a 2-d model for which accurate numerical results are available [|l^. This will allow us to 



compare the results obtained by incorporating subcritical bubbles into the calculation of the 
decay barrier with the results from the numerical simulations. 

The 2-d scalar potential V{(j)) is given in eq. (|I|). Following the rescaling of Ref. ||1Q|| , 
the potential can be written in terms of one dimensionless parameter A = rn^h/ g'^, 

V{<P) = \<P' - + (13) 

This double-well potential is degenerate when A = 1/3, and the second minimum is lower 
than the first when A < 1/3. 

As argued before, we find the new coarse-grained potential Kg (or, equivalently, Acg) by 
constraining it to agree with the thin wall limit. Simple algebra from eqs. (^ and (|13]) yields, 
to first order in the deviation from degeneracy, 

^c. = A - If (14) 

~ 2 

where T^c = ^^sc is the dimensionless free-energy density in subcritical bubbles. The new 
potential Vcg is then used to find the bounce solution and the free energy of the critical 
bubble. 
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The calculation of J^sc in two dimensions is fairly straightforward. Close to the thin 

wall limit {i.e., Go^+ ~ Gj, ,0 = G), one can analytically solve the equilibrium Boltzmann 

equation for the density distribution function, obtaining f{R,(f)A,T) = (1 — 2'j)Wt{R, 4>a), 
where {v = \v\) 



WTiR,<f>A) = ^^exp 




-- + RT{a/v)+ ^ ' 




(15) 



and we wrote the free energy of a given subcritical configuration as = a + I3R?, with 
a = 7r(j)\/2, and f3 = a (^^ — |</)^ + • The fraction of the volume occupied by subcritical 

bubbles is then, 

/ 1 T3 D \ -^T(0max; -Rmin; -Rmax) , 

Tl'Pmax) -ttmiri) -ttmaxj — TTj , t / i d d vT ' [^^) 

[I + ^-iT(,<Pmax5 -n-min, -n-maxjj 

where. It = I^^^ /r™'"' i^r'^WTdrdcj) . The radial integration can be done analytically, 
although the result is not particularly illuminating. The integral over amplitudes must 
be done numerically. We then substitute f{R,(f)A,T) and 7 into eq. (p!2[) to finally find 
J^sc(A,T, A/v,a/v). 

In order to illustrate the effect of subcritical bubbles on the nucleation barrier, in Figure 
1 we compare the value for the barrier obtained with and without the corrections, as a 
function of A, for constant values of the temperature. The temperatures are chosen to be 
within the range used in the 2-d simulation. The constant A was fixed aX A = 0.02, consistent 
with the measurements of Ref . |10| . Notice that the presence of subcritical bubbles greatly 
decreases the barrier as the potential approaches degeneracy (A —>■ 1/3). However, for small 
temperatures T < lOm^/g'^, the correction becomes negligible. 

In Figure 2 we show that the calculation of the nucleation barrier including the effects 
of subcritical bubbles is consistent with data from lattice simulations, whereas the standard 
calculation overestimates the barrier by a large margin. In fact, the inclusion of subcritical 
bubbles provides a reasonable explanation for the anomalously high nucleation rates observed 
in the simulations close to degeneracy. The error bars are from the numerical measurements 
of the barrier; for larger values of A, higher temperatures had to be used to attain nucleation, 
increasing the error in the barrier measurements. However, we note that even with the large 
error bars, the data is inconsistent with the theoretical predictions for the barriers, while the 
corrected barrier values fall within the error bars for a wide range of parameters. We note 



that data from 1-d simulations also show the same behavior as the data in Figure 2 |jTO|. 
Simulations in 3-d are in progress, and will enable us to test this method in more detail. 

Finally, we note that the inclusion of nonperturbative corrections through the definition 
of an effective, "coarse-grained", coupling may have several consequences not only to the 
nucleation rate of first order transitions, but also to their dynamics. Clearly, once we have a 
corrected potential, quantities such as the critical temperature, the amount of super-cooling, 
the bubble-wall velocities, and the completion time for the transition will change. This 
opens up several possible applications of this method, from laboratory studies of nucleation 
to cosmological phase transitions. 
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List of Figures 

Figure 1. Comparison of the decay barrier as a function of A with and without the in- 
clusion of subcritical bubbles, at fixed temperatures. 

Figure 2. Comparison between numerical data and theoretical predictions for the decay 
barrier with and without the inclusion of subcritical bubbles. 
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